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Abstract 

When the cross section of J/'ij) production is considered varying with the energy of the 
nucleon-nucleon interaction the production of J/ip in pA and A A collisions has been stud¬ 
ied using FRITIOF Model. The calculation shows that the cross section oi J/^j) production 
“per nucleon-nucleon collision” decreases with increasing mass number and centrality as a con¬ 
sequence of continuous energy loss of the projectile nucleons to the target nucleons in their 
successive binary nucleon-nucleon collisions. We have compared our model predictions with the 
experimental data of J/i/: production. 
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Suppression oi Jftp production in high energy heavy ion collisions was proposed as an effective 
signature of QGP formation ten years ago [Q . The ensuing experimental data confirmed a significant 
suppression of J/^/> production in both pA and AA collisions However, alternative explanations 
of the J/ip suppression exist based on the absorption of J/ip in nuclear matter |^. The overall set 
of extensive data collected and analysed by NA38 seems to support the absorption mechanism in 
P+Btarget, 0+^target and S+U collisions. Nevertheless a rather large absorption cross section 
6.2 mb) has to be used in order to fit the experimental data, about three times larger than the total 
J/'0-N cross section from the EMC Collaboration Such a difference was already noted long time 
ago in . Recent calculations based on the colour octet model show that this phenomenological cross 
section could be understood by the absorption of pre-resonance states {ccg) in nuclear matter |^. 
Other sources oi Jjip suppression in nuclear collisions, such as the interaction of J/ip particles with 
the produced mesons (called comover) j^, gluon shadowing in nuclei, intrinsic charm component, 
energy degradation of produced cc pair, etc. |^, have been introduced besides the absorption due 
to the J/ip-N interaction to explain the data. 

We in this paper propose a simple model to investigate J/ip suppression in pA and A A collisions 
focusing on the decrease of the cross section oi J/ip production with increasing mass number and 
centrality due to the continuous energy loss of the projectile nucleons to the target nucleons in their 
successive binary nucleon-nucleon collisions, not on its later absorption in nuclear matter. From 
the calculations of this model, we conclude that the absorption oi J/ip particles in nuclear matter 
is only part of the sources of J/ip suppression seen in pA and AA collisions so far. 

Since the probability for a nucleon-nucleon collision leading to J/ip production is very small it is 
generally accepted that multiple J/ip production processes in pA and AA collisions can be neglected. 
The studies oi J/ip suppression based on the absorption mechanism actually assume that the prob¬ 
ability oi J/ip production “per nucleon-nucleon collision” is the same independent of the masses of 
the colliding nuclei at a given energy. However, each binary nucleon-nucleon collision experienced 
by a projectile nucleon on its way out of the target in pA and AA collisions may not be the same 
if the projectile nucleon loses a fraction of its energy in each binary collision to the target nucleon 
such that the cms energy of each binary collision of this incoming nucleon with a target nucleon is 
different. For the cross sections of the hard QCD parton-parton scatterings, which is responsible 
for J/ip production, increase with increasing energy it is conceivable that the probability of J/ip 
production would also depend on the energy of a nucleon-nucleon collision in a similar way. It is 
confirmed both theoretically and experimentally that J/ip photoproduction cross section exhibits a 
strong threshold behaviour at the low energy region and then increases with energy at the higher 
energy region [||. 

In a participant-spectator model of nucleus-nucleus collisions, like FRITIOF, each projectile nucleon 
may collide several times with the target nucleons. If momentum transfers are assumed to take place 
in each binary nucleon-nucleon collision, then the cms energy of these binary collisions will decrease 
with time, thereby the probability to produce J/ip in each binary collision going down. 

We in this paper have calculated how the cross section oi J/ip production varies with increasing mass 
number of colliding nuclei and centrality using FRITIOF model. The calculations show that J/ip 
production is suppressed in pA and AA collisions in comparison with the nucleon-nucleon collision, 
and the larger the mass number of colliding nuclei and the centrality, the greater the suppression, as 
a result that probability of finding a QCD hard process in a binary nucleon-nucleon decreases with 
increasing mass number of colliding nuclei and centrality. Taking into account the J/ip absorption 
by J/ip-N interaction our results are in good agreement with experimental data with the exception 
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of the latest data of Pb+Pb collisions at 158 AGeV/c, which show a further J/ij; suppression [ p^ . 
The absorption cross section aabs needed to explain the data from p+Btar-get, 0+^target and S+U 
collisions is about 1.4 mb in this paper, which is consistent with the experiments (EMC Collaboration 
0 gives a total J/ili-H cross section 2.2 ± 0.7 mb and the J/'i/l’-N quasi-elastic cross section is given 
in as 0.79 ± 0.012 mb). Furthermore, our results also imply that some new mechanism of Jhb 
suppression seems to be needed particularly for understanding Pb-|-Pb data. The authors in Q 
have attributed the further J/ip suppression in Pb-fPb data to the formation of a QGP state. 

As we have mentioned before the cross section of J/ip production will be different in each binary 
nucleon-nucleon collision in pA and AA collisions if the energy loss of the projectile nucleons to the 
target nucleons in the successive binary collisions is taken into account. Assume that < > is 

the mean cross section for the production of a J/ip particle in a binary nucleon-nucleon collision 
(here the average is done over all the binary collisions at an impact parameter h), then the total 
probability for producing & J/ip particle in the collisions of A -|- B at an impact parameter b is the 
sum [Q 

A.B / ^ ^ \ 

= n )™>ni-(1) 

where T(b) is the thickness function. Because T(h) < > is a very small quantity (< 

10“"*’ fm^ 1^, r(6) ~ 10“^fm“^ for central S-l-U collisions, for instance), the summation given 
by Eq. i) is dominated by the first term with n=l. The terms with n>l represent multiple J/ip 
production processes and shadowing corrections, which are very small and can be neglected. Then 
the probability for J/ip production in A-l-B collisions can be approximated to be 

Pj/%=AB[T{h)<a%l>]. (2) 

Therefore, the cross section of J/ip production corresponding to a centrality bin, A5, is given by 
the following formula 

= AB< aj/l >Ab [ mdb, (3) 

J A6 

where < >Ab is the mean cross section of J/ip production “per nucleon-nucleon collision” 

within the centrality bin A6. We know that QCD hard scatterings (the gluon fusion and quark- 
antiquark annihilation) between partons are the main source of J/ip production. Let Ph be the 
probability to have a hard scattering in a nucleon-nucleon collision and the probability to 

produce a J/ip from the hard scattering, then we can write out < > to be 

< >=< ^tPu^^Ph >= ^tPI''^ < Ph >, (4) 

where ut is the total cross section of a nucleon-nucleon collision. We have assumed that arPj^^^ is 
approximately a constant in the energy span that we are concerning, so the product is the same for 
all the binary collisions. Combining Eq.(^) with Eq.(||) and replacing < Ph > hy the ratio 

{nbinib) is the number of binary collisions in an A-l-B collision and rihini^) number of binary 
collisions with a hard scattering. Both of them are the function of the impact parameter b) we 
finally obtain 

= ABarP^i'* < >Ab [ T{b)db (5) 

BhinXp) J Ab 


and for the minimum bias events we have 






We see from Eq.(||) that the dependence of the quantity on the masses of colliding nuclei or 
centrality is solely determined by how the mean probability of having a hard scattering in a binary 
nucleon-nucleon collision varies with the masses of colliding nuclei or centrality. Before calculating 
in pA and AA collisions we will give a brief introduction of FRITIOF dynamics focusing on 
how a hard parton-parton scattering is distinguished from a soft one. 

FRITIOF is a string model based on the concepts of the Lund String Model |^, which started 
from the modeling of inelastic hadron-hadron collisions and it has been successful in describing 
many experimental data from the low energies at the ISR-regime all the way to the top SPS ener¬ 
gies [0 0 - This has been achieved by the introduction of a particular longitudinal momentum 
transfer scenario, gluon bremsstrahlung radiation (The Dipole Cascade Model, PC M |1^ , and the 
Soft Radiation Model, SRM [Q, — this is implemented by the use of ARIADNF ^|) as well as 
hard parton scattering (Rutherford Parton Scattering, RPS — this is implemented by the PYTHIA 
routines [^). In FRITIOF, during the collision two hadrons are excited due to longitudinal mo¬ 
mentum transfers and/or a RPS. It is further assumed that there is no net color exchange between 
the hadrons. The highly excited states will emit bremsstrahlung gluons according to the SRM. They 
are afterwards treated as excitations or the Lund Strings and the string states are allowed to decay 
into final state hadrons according to the Lund prescription as implemented by JFTSFT . 

In the FRITIOF model a hadron is assumed to behave like a massless relativistic string (MRS) 
corresponding to a confined color force field of a vortex line character embedded in a type II color 
superconducting vacuum. A hadron-hadron collision is pictured as the multi-scatterings of the 
partons inside the two colliding hadrons. This includes both the hard and the soft components 
depending on the four-momentum transfers Q^, or equivalently the transverse momentum transfers 
involved. The soft part is described by a simple phenomenological model. The hard scatterings can 
however be calculated from perturbative QCD, and correspond to the Rutherford parton-parton 
scattering (RPS). The divergence problem in RPS is handled by introducing the Sudakov factor. 

There will be color separation in the model, i.e. there will for each hadron be a color 3 (a “diquark”) 
continuing forward along the beam direction and a valence quark, a color 3, moving in the opposite 
direction due to the longitudinal momentum transfer. This will lead to bremsstrahlung of a dipole 
character. 

A procedure therefore is adopted in FRITIOF that compares the “hardness” of the Rutherford 
partons to that of the bremsstrahlung gluons. The RPS is accepted only if it is harder than the 
associated radiation. If the RPS is “drowned”, which is to say that it is softer than the radiation, 
then the RPS is not acceptable and the collision proceeds as a purely soft collision. With this 
prescription the RPS spectrum is suppressed smoothly at small to medium transverse momentum 
region. 

For the hadron-nucleus and nucleus-nucleus collisions, the process has in the FRITIOF model been 
treated as a set of incoherent collisions on the nucleons. Thus a nucleon from the projectile interacts 
independently with the encountered target nucleons as it passes through the nucleus. The probability 
distribution for the number of inelastic collisions ly is taken from geometric calculations. Fach of 
the sub-collisions is treated in the same way as an ordinary hadron-hadron collision, although 
the momentum transfers will again be additive and every encounter will make the projectile more 
excited. If it interacts with z/ nucleons in the target, -I- 1 excited string states will be formed as 
a result. These string states will then independently emit associated bremsstrahlung radiation and 
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then fragment into hadrons in the same way as individual strings. This picture is supported by 
the fact that the global features of heavy ion collisions are satisfactorily explained by the collision 
geometry together with the independent hadron-hadron collisions. 

Using FRITIOF it is straightforward to calculate the number of binary nucleon- nucleon collisions, 
nbin{b), and the number of the binary collisions with a hard scattering, in pA and AA 

collisions at a given impact parameter &, so that the mean probability to have a hard scattering in a 
binary nucleon-nucleon collision, < Ph > = , can be obtained. Since we are mainly interested 

in J/Ip production in pA and AA collisions relative to that in the pp collision we do not need to 
know how a J/ip is actually formed from the hard scatterings in order to investigate the dependence 
of J/ip production cross sections on mass number and centrality. 

We have calculated the quantity for various pA and AA collisions (and < >Ab for 

different centrality bins in S-l-U and Pb-fPb collisions) at P/ab=200 GeV/c using FRITIOF. After 

, ^AB 

determining (TtP/^ ^ by the data of the pp collision we plot our results of as a 

function of Aproj. ■ Btarg. in Figure 1 for the minimum bias events. For the cross section of J/ip 

^AB 

production in different centrality bins we plot the results of as a function of Np-Nt in Figure 

2, where Np is the number of participants from the projectile and Nt the number of participants 
from the target, since the Drell-Yan cross section in a given centrality bin is found in experiments 
proportional to an effective Aproj. ■ Btarg. (= Aproj. ■ Btarg. T(b)db). In the same way, a constant 
has to be determined by the corresponding data of the pp collision. The impact parameter bins 
are taken to be the same as those extracted by NA38 and NA50 |^. We decided not to use 
the absorption length L to be the longitudinal axis as used by NA50 because L calculated from 
the geometry model is not sensitive to the change of impact parameter for very central Pb-|-Pb 

^AB 

collisions. The results of our calculations show that the decrease of quantity, Bpp^ - ^ - or 

P'P' AAp-pQj_ ' iJtarg . 

^AB 

Bpp-A^ (the cross section oi J/ip production “per nucleon-nucleon collision”), is due to the fact 
that the probability of the QCD hard scattering per binary nucleon-nucleon collision decreases with 
the increasing mass number and centrality. 

When the absorption of J/ip by J/ip-H interaction is also taken into account, i.e. the previous 
results are multipied by exp(—pLaabs) with p=0.14 n/fm^, aabs=^A mb and L taken to be the same 
as those in |^, our model reproduces the data of J/ip suppression with the exception of the latest 
data from Pb-|-Pb collisions at 158 GeV/c, which clearly show a further suppression. 

One possibility which will bring about a further suppression of J/ip production is that the 

probability to produce a J/ip from a hard process, drops down suddenly under certain conditions. 
This is equivalent to say that the cc produced from the hard processes can not form a bound state. 
However at the moment our simple model can not estimate when this would happen. 

However, there are still other possible mechanisms of J/ip suppression, which are not included in 
our simple model. The xp— dependence of J/ip suppression is not investigated yet. Therefore, it 
is hard to make any conclusion now whether this further suppression in Pb-fPb collision is due to 
QGP formation. 

It is known that there is no unique criterion to distinguish a hard process from a soft one in a nucleon- 
nucleon collision. Usually a qTmin is introduced to be the minimum transverse momentum of the 
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produced partons from a hard process. A dynamic criterion is applied in FRITIOF to chose a hard 
process by comparing the hardness of a RPS parton with the hardness of the bremsstrahlung gluons 
as mentioned before. However, the cross section of J/ip production should not depend on which 
criterion is actually used in the calculation. We have thus calculated all the results of this paper 
using the conventional qTmin criterion in Pythia {qTmin= IGeV/c), just to check if our conclusion 
relies on the specific criterion in FRITIOF. The calculations show that the results in these two cases 
are in agreement with each other. 

We have also checked if the quantity is energy-independent as we have assumed. A 

parametrization form of the J/tp cross section is given as | p4[ 

(7) 

VS 


where \fs stands for the cms energy per nucleon. Therefore, if in Ed-® is energy- 

independent then we should have a ratio 


(1 





'^bin 


= C = 


constant. 


( 8 ) 


We have calculated 

'^bin 

and the ratio at various energies from P;jj{,=60 GeV/c to Piab= 450 GeV/c for pp collisions {ni,in= 
1 for a pp collision) and the results are listed in Tab.l, which show that this ratio in this energy 

. . . rt'*. 

region is not sensitive to the change of energy in comparison with . But a threshold behaviour 
may exist at lower energies, which can be seen from the value at Piah=60 GeV/c. 
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Figure Captions 


Figurel: J/'0 cross sections divided by Aproj. • Btarg. as a function of ^roj. • Btarg.- Our results are 
compared with the experimental data from NA51, NA38 and NA50 [|l0[. 


Figure2: J/ijj cross sections divided by Drell-Yan cross sections as a function of Np ■ Nt- Our results 
are compared with the experimental data from NA38 and NASO 0- 
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Table 1. The energy dependence of the probability to 
have a hard scattering per binary collision 
and the assumed constant C in Eq.(||) 


Piab GeV/c 


bin 

G 

60 

1.17 

X 

10-^ 

1.84 

100 

2.82 

X 

Tir^ 

2.03 

150 

4.76 

X 

Tir^ 

2.14 

200 

6.40 

X 

Tir^ 

2.22 

300 

9.17 

X 

Tir^ 

2.27 

450 

12.2 

X 

Tir^ 

2.30 















^(J/^)/(AprojectiieBtarget)(nb/nucleon) 


* p (450GeV/c ) - A (A = P, d, Cu, W) (NA51, NA38) 



A(projectile)B(target) 





